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ABSTRACT 

High contrast coronagraphic imaging of the immediate surrounding of stars requires exquisite 
control of low-order wavefront aberrations, such as tip-tilt (pointing) and focus. We propose an 
accurate, efficient and easy to implement technique to measure such aberrations in coronagraphs 
which use a focal plane mask to block starlight. The Coronagraphic Low Order Wavefront Sen- 
sor (CLOWFS) produces a defocused image of a reflective focal plane ring to measure low order 
aberrations. Even for small levels of wavefront aberration, the proposed scheme produces large 
intensity signals which can be easily measured, and therefore does not require highly accurate 
calibration of either the detector or optical elements. The CLOWFS achieves nearly optimal 
sensitivity and is immune from non-common path errors. This technique is especially well suited 
for high performance low inner working angle (IWA) coronagraphs. On phase-induced amplitude 
apodization (PIAA) type coronagraphs, it can unambiguously recover aberrations which origi- 
nate from either side of the beam shaping introduced by the PIAA optics. We show that the 
proposed CLOWFS can measure sub-milliarcsecond telescope pointing errors several orders of 
magnitude faster than would be possible in the coronagraphic science focal plane alone, and can 
also accurately calibrate residual coronagraphic leaks due to residual low order aberrations. We 
have demonstrated sa 10~ 3 A/D pointing stability in a laboratory demonstration of the CLOWFS 
on a PIAA type coronagraph. 

Subject headings: instrumentation: adaptive optics — techniques: high angular resolution 



1. Introduction 

High performance coronagraphs with small in- 
ner working angle (IWA) are unavoidably very 
sensitive to smal l pointing errors and other low or- 
der aberrations (TLlovd h Sivaramakrishnanl 2005; 
Shaklan fc Greenl 120051: " ISivaramakrishnan et al 



20051 : iBelikov et all I2006t iGuvon et alj |2o63T 

This property is due to the fact that the wavefront 
of a source at small angular distance (typically be- 



tween 1 and 2 X/D) from the optical axis is "sim- 
ilar" (in the linear algebra sense of the term) to 
an on-axis wavefront with a small (<C X/ D) point- 
ing error. If the coronagraph must "transmit" the 
former, it will also transmit a significant part of 
the latter, and therefore be extremely sensitive to 
pointing errors. This behavior is indeed verified 
by perfor mance comparison between coronagraph 
concepts ( Guvon et al.ll2006h . 
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Fig. 1. — Optical layout of a coronagraphic low order wavefront sensor system, shown here with a PIAA 
coronagraph. See text for details. 



For high performance coronagraphs, such as the 
Phase-Induce d Ampl i tude Apodization (PIAA) 



coronagraph ( Guvonl 2003 : Guvon et al. 2005 ) 



used as example in this paper, the coronagraph 
should ideally be designed to balance IWA against 
stellar angular diameter, which sets a fundamental 
limit on the achievable coronagraph performance. 
Such coronagraphs are therefore pushed to be- 
come sensitive to pointing errors corresponding to 
the angular size of nearby stars, roughly 1 mil- 
liarcsecond (mas), and are also highly sensitive to 
other low order wavefront errors such as focus and 
astigmatism. Milliarcsecond-level pointing error 
can increase stellar leakage in the coronagraph to 
the point where a planet would be lost in photon 
noise. Even smaller errors can create, if not indc- 
pendantly measured, a signal which is very similar 
to a planet's image. 

A robust and accurate measurement of low or- 
der aberrations (especially tip-tilt errors, which 
are easily generated by telescope pointing errors 
and vibrations) is therefore essential for high con- 
trast coronagraphic observations at small angu- 
lar separation. The science focal plane after the 
coronagraph is unfortunately " blind" to small lev- 



els of low order aberrations, which can only be 
seen when already too large to maintain high con- 
trast in the coronagraphic science image. A bet- 
ter option is to monitor pointing errors by us- 
ing starlight which would otherwise be rejected 
by the coronagraph. This scheme was successfully 
implemented on the LYOT project coronag raph 



(Op penheimer et alJl2004HDigbv et alJl2006f l. Al- 



ternatively, the measurement could be performed 
independently from the coronagraph optical train 
(for example, the wavefront sensor in the adap- 
tive optics system upstream of the coronagraph). 
We propose in this paper an improved solution 
to obtain accurate measurement of several low- 
order aberrations including pointing: the coron- 
agraphic low-order wavefront sensor (CLOWFS). 
The wavefront control requirements for a PIAA 
coronagraph are first clearly defined in 321 The 
CLOWFS principle is presented in $3] Wavefront 
reconstruction algorithms and CLOWFS sensitiv- 
ity are discussed in 21 The aberration sensitivity 
of a PIAA coronagraph equipped with a CLOWFS 
is discussed in $5l and the results of a laboratory 
demonstration on a PIAA Coronagraph system 
are shown in S}6l 
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Low order wavefront control require- 
ments for PIAA coronagraphs 
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Fig. 2. — Simulated science focal plane image in a 
PIAAC system with a 0.005 A/D pointing offset. 
The PIAA design adopted for this simulation has 
a 1.9 A/D 1WA. 

Wavefront aberrations can be produced either 
before the PIAA optics (for example a bending of 
the telescope primary mirror, or a telescope point- 
ing error) or between the PIAA optics and the fo- 
cal plane mask. In a space-based telescope free of 
atmospheric turbulence, the strongest sources of 
aberration are likely to be telescope pointing er- 
rors and low-order aberrations due to structural 
deformations of the optical telescope assembly. 
The PIAA coronagraph is especially sensitive to 
such low order aberration s if they are introdu ced 
prior to the PIAA optics (Beli kov et al. 2006f ). as 



they will then scatter light in the most scientifi- 
cally precious area of the science focal plane: the 
inner part of the field. As illustrated in Figure 
[21 a very small pointing error (0.5% of A/D) can 
be sufficient to create an artefact as bright as an 
Earth- like planet. Extreme sensitivity to pointing 
is unavoidable in small IWA coronagraphs, and, 
in the case of the PIAA, is due to the fact that a 
"remapped" tip/tilt scatters light outside the focal 
plane mask. Coronagraphs with larger IWA and 



better tolerance to pointing errors exist, and even 
within the PIAA "family" of coronagraph, sensi- 
tivity to pointing errors can be balanced against 
IWA by changing the size of the focal plane mask 
and the pupil apodization profile. 

The PIAA coronagraph is much less affected by 
aberrations after the beam shaping optics. For ex- 
ample, a small amount of tip/tilt in the post-PIAA 
optics will simply move the PSF on the focal plane 
mask (which is most likely slightly oversized) with- 
out introducing scattering at larger separations. 
In addition, low order aberrations introduced af- 
ter the PIAA optics will likely be much smaller 
in amplitude, thanks to the small size of optical 
elements between the PIAA optics and the focal 
plane mask. 

For successful detection of a faint source (as- 
sumed here to be an Earth- like planet at 10~ 10 
contrast), low-order aberrations must simultane- 
ously: 

• Be small enough to avoid loosing the planet 
signal within the scattered starlight's pho- 
ton noise. For the PIAA coronagraph de- 
sign used in this work, a pre-PIAA 0.005 
A/D pointing error is sufficient to scatter 
6.9 x 10~ 10 of the starlight into the science 
field, and this scattered light produces two 
wide "arcs" on either side of the optical axis, 
at a 10 -10 contrast level. Although a PIAA 
coronagraph could be designed with reduced 
sensitivity to pointing errors (but with larger 
IWA), we assume here that the maximum 
allowable pointing error for a space PIAA 
coronagraph mission aimed at direct imag- 
ing of "Earth-like" planets is 0.005 A/D (cor- 
responding to starlight leak peaking at 10~ 10 
contrast). On a 1.4-m diameter telescope in 
the visible, this corresponds to 0.4 mas, a 
value which is similar to the angular radius 
of a "typical" target (a main-sequence star 
at 10 pc). On larger telescopes, the angular 
radius of the star (0.5 mas for a Sun-like star 
at 10 pc) will drive the coronagraph design, 
which will therefore also end up with a s=s 0.5 
mas RMS pointing error requirement. 

• Be stable, or calibrated, to a fraction of the 
planet's expected flux. We assume in this 
paper that coronagraphic leaks due to low 
order aberrations must be calibrated to 10% 
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of the expected planet's contribution. This 
second requirement is therefore much more 
severe. Light scattered in the science focal 
plane scales approximately as the square of 
the aberration amplitude: pointing must be 
stable (or calibrated) to « 0.0016A/D (0.13 
mas on a 1.4-m telescope in the visible) for a 
planet at 1CP 10 contrast. Fortunately, high 
accuracy measurement of low order aberra- 
tions by the CLOWFS scheme proposed in 
this paper can be used to reliably model the 
stellar leakage in the coronagraphic science 
focal plane. This model can then be numer- 
ically subtracted from the science image to 
reveal much fainter underlying sources down 
to the photon noise and detector noise lim- 
its. Pointing errors larger than 0.0016 X/D 
(but smaller than 0.005 X/D) are acceptable, 
as long as they are measured to 0.0016 X/D 
accuracy. 

3. Coronagraphic Low Order Wavefront 
Sensor: Principle 

3.1. Optical Layout 

A simplified optical layout for a CLOWFS sys- 
tem for a PIAA coronagraph is shown in Figure 
[TJ CLOWFS light is extracted by the focal plane 
mask located after the PIAA optics. In the PIAA 
coronagraph design, the role of this focal plane 
is to selectively remove starlight, while transmit- 
ting the science field. The mask is therefore illu- 
minated by a large number of photons, which are 
freely available, and, if properly used, allow highly 
sensitive measurement of low order aberrations. 

The central part of the focal plane mask used in 
the CLOWFS design is opaque: only a reflective 
annulus around this central part sends light to the 
CLOWFS optics. As shown in figure [T] we denote 
r\ the radius of the inner opaque zone and r-i the 
radius of the focal plane mask, which is fixed at 
the IWA of the coronagraph. For the baseline con- 
figuration adopted in this paper r\/r% = 0.4, with 
r 2 w 1.8X/D on the sky: the central 40% of the 
focal plane mask is opaque. Since the pupil be- 
fore the focal plane mask is strongly apodized by 
the PIAA optics, only a small fraction of the total 
light is reflected by the reflective annulus covering 
the r\ < r < r 2 range (see Figured]). 

The CLOWFS detector acquires a defocused 



image of this annulus, at defocus distance 20rad — 
3.2 waves (in this paper, the defocus value is given 
as Peak-to- Valley in the post-PIAA pupil plane). 
While this defocus value may seem large, most 
of the light is in the central part of the apodized 
pupil, and the "effective" defocus is about half to 
a third of the optical defocus listed in this work. 
The CLOWFS pupil is oversized by a factor 2 
(/3 = 2a) to include most of the light reflected 
by the CLOWFS focal plane mask. 
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Fig. 3. — Fraction of the total starlight collected 
by the telescope sent to the CLOWFS detector as 
a function of r% /r<i . 

3.2. Importance of defocusing the CLOWFS 

The CLOWFS defocus (the CLOWFS detector 
is not conjugated to the focal plane mask) is nec- 
essary to measure focus, but is not necessary if 
only tip and tilt are measured. In a CLOWFS 
configuration where the detector is conjugated to 
the focal plane mask, a single CLOWFS image 
can be used to recover the amplitude of focus in 
the incoming beam, but not its sign (in an opti- 
cal system free of aberrations, focal plane images 
acquired inside and outside focus are identical). 
The CLOWFS defocus is therefore introduced to 
remove this sign ambiguity. A mathematical proof 
of the quadratic (rather than linear) response of a 
"in- focus" CLOWFS to pupil focus aberration is 
provided in 34. II 

3.3. Why a dual zone focal plane mask ? 

While in an ideal system (where photon noise 
would be the only source of noise) the central 
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opaque zone would reduce the CLOWFS perfor- 
mance, it is in practice extremely advantageous 
for two reasons, which are now described. 

3.3.1. Relative signal amplification 

Without the opaque zone, the CLOWFS de- 
tector would be illuminated by a large number of 
photons and the signal to be extracted would be 
a tiny relative change in intensity, which would be 
challenging to measure in practice. To detect this 
signal, the detector calibration would need to be 
very accurate and detector saturation would be 
a serious issue. Masking the central part of the 
PSF greatly reduces the total amount of light in 
the CLOWFS but only slightly reduces the ampli- 
tude of the signal produced by low order aberra- 
tions. For example, a small pointing error in the 
post-PIAA pupil creates maximal signal where the 
PSF surface brightness slope is the greatest. The 
very center of the PSF, although it contains a lot 
of flux, contains very little signal. Thanks to the 
central dark zone, microscopic changes in the coro- 
nagraphic wavefront produce macroscopic changes 
in the CLOWFS image, as shown in Figured! 

3.3.2. Self -referencing 

The central dark zone of the focal plane mask 
allows the CLOWFS to be insensitive to small mo- 
tions of the CLOWFS elements (detector and off- 
axis parabola or lens to re-image the focal plane 
mask) . 

If the focal plane mask were fully reflective, a 
tip-tilt error in the post-PIAA pupil (which should 
be corrected) and a tip-tilt error in the CLOWFS 
optics (which should not be corrected) would look 
identical. This is due to the fact that the outer 
edge of the focal plane mask is dark, and the 
CLOWFS would therefore measure tip-tilt as a 
translation of the image on the CLOWFS detec- 
tor. This sensitivity to small amount of tip-tilt in 
the CLOWFS optics would therefore require the 
CLOWFS optical elements positions to be accu- 
rately known and stable. With the dark zone in 
the center of the focal plane mask, however, the 
same tip-tilt is measured as a macroscopic change 
in the CLOWFS image shape (see Figure 0]) in- 
stead of a small translation. Thanks to the dark 
zone, the CLOWFS tip-tilt measurement is accu- 
rately referenced to the focal plane mask. 
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Fig. 4.— Images obtained by the CLOWFS for 
T\jr<x — (left) and ri/ra = 0.4 (right) for a 
defocus distance of 20rad. With n/r2 — 0.4, 
the images are significantly fainter but aberrations 
are more easily seen. In each of the 8 images 
shown, the inside and outside focus images of the 
CLOWFS mask are side-by-side. In the CLOWFS 
system we propose, only one of these two images 
would need to be acquired. 



3.4. Sensitivity to focal plane mask man- 
ufacturing errors 

The geometry of the focal plane mask is driven 
by both the coronagraph architecture and the 
CLOWFS. The transmission of the mask must 
first satisfy the coronagraph's requirements. Im- 
plementing a CLOWFS requires: 

• A tilted focal plane mask to redirect stellar 
light to the CLOWFS imaging optics. If the 
tilt angle is large, or if the coronagraph re- 
quires a very circular focal plane mask, the 
focal plane mask could be made elliptical. 

• A reflective zone on the focal plane mask. 
This requirement should have no negative 
impact on the coronagraph's performance, 
as a reflective coating may be deposited on 
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the illuminated side of the focal plane mask 
without affecting its transmission. 

The focal plane mask reflectivity defines the 
CLOWFS performance. An "ideal" focal plane 
mask for CLOWFS, as described in this paper, 
may be very challenging to manufacture: 

• ft is very difficult to make the central part 
of the mask truly " black" , and some of the 
light in the central part of the PSF will be 
reflected into the CLOWFS imaging optics. 

• The reflective annulus may not be uniformly 
reflective and could also introduce wavefront 
errors (the reflective surface may not be very 
flat). 

• Coatings, wether optimized to be black (cen- 
tral part of the mask) or reflective (annulus), 
are somewhat chromatic: the reflectivity 
map of the focal plane mask is wavelength- 
dependent. 

Even with these imperfections, the CLOWFS will 
still produce a strong response for the aberrations 
to be measured (tip, tilt, focus and their remapped 
equivalents). It may however be difficult to pre- 
dict what there responses will be if the focal plane 
mask is poorly calibrated. As detailed in fJU we 
therefore propose to first measure these linear re- 
sponses by introducing aberrations in the coro- 
nagraph optics, and then use these responses to 
decompose the CLOWS image in a linear sum of 
aberrations. This step requires no additional hard- 
ware, provided that there are actuators to correct 
for the aberrations measured by the CLOWFS. 
This is very similar to what is commonly done on 
adaptive optics systems, where a "response ma- 
trix" is first acquired by measuring the WFS sig- 
nal when each actuator of the deformable mirror 
is moved. 

4. Wavefront reconstruction algorithm 
and sensitivity 

4.1. Linearity 

We show in this section that for small wavefront 
aberrations (<< 1 rad), the CLOWFS image is a 
linear function of the wavefront aberration modes 
to be measured. This convenient property is due 
to: 



• The fact that the CLOWFS is not operat- 
ing on a "dark" fringe: when the incoming 
wavefront is perfect, the CLOWFS images 
already contain some light. 

• The non-orthogonality between the aberration- 
free CLOWFS complex amplitude distribu- 
tion in the CLOWFS detector array and the 
change introduced on this complex ampli- 
tude by the aberration modes to be mea- 
sured. 

We denote u and x the 2-D coordinate respectively 
in the post-PIAA pupil plane and in the CLOWFS 
detector plane. We denote Aq(x) the 2-D com- 
plex amplitude obtained by the CLOWFS for a 
perfect wavefront W(u) — 0. The complex am- 
plitude A(x) obtained in the CLOWFS detector 
plane is a linear function of the pupil plane com- 
plex amplitude (virtually all optical systems are 
linear in complex amplitude). Since by definition, 
A(x) — Aq(x) for W(u) = 0, we can therefore 
write 

A{x) = Aq(x) + MW(u). (1) 

where M is a linear operator. The correspond- 
ing light intensity I(x) = \A(x)\ 2 in the CLOWFS 
detector plane is 



I(x) = I (x) + 2Re A (x)MW(u) + \MW(u)\ 2 

(2) 

which is linear function of W(u) as long as : 



\MW{u)\ z < Re A {x)MW(u) 



(3) 



For small wavefront aberrations, |MW(u)| <C 
Ao(x)|, which helps to satisfy condition [3] If 
there is an aberration mode for which Aq(x) and 



= 



MW(u) are orthogonal (Re \a (x)MW(v) 

for all values of a;), condition[3]will not be satisfied 
even for small aberration levels, and the CLOWFS 
image I(x) will be a quadratic function of the 
wavefront aberration (I(x) — Iq(x) + \MW(u)\ 2 ). 

Although such a situation is " unlikely" because 
orthogonality would need to occur simultaneously 
on each pixel of the CLOWFS detector, it does oc- 
cur for the focus aberration mode if the CLOWFS 
detector is conjugated to the focal plane mask (no 
defocus in the CLOWFS re- imaging) . In this spe- 
cial case, since the CLOWFS re-images the fo- 
cal plane reflective annulus without defocus, we 
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can look at the orthogonality between A (x) and 
MWF(u) directly on the focal plane mask. The 
complex amplitude Aq (x) is the Fourier transform 
of the pupil plane complex amplitude and is purely 
real (no imaginary part) over the reflective annu- 
lus of the focal plane mask when no wavefront 
aberration is present. MW(u) is the change in 
the complex amplitude over the focal plane annu- 
lus introduced by a focus aberration in the pupil 
plane, and is therefore the Fourier transform of 

(l-e*M)/(r)«tp(r)/(r) (4) 

where r is the radial coordinate in the pupil plane, 
f(r) is the amplitude profile in the exit pupil of 
the PIAA coronagraph (/(r) decreases with r), 
and p(r) is the phase in the exit pupil of the 
PIAA coronagraph (p(r) is the Focus phase term 
remapped radially by the PIAA optics). The ap- 
proximation in Equation U is valid because we are 
considering small aberrations (p(r) << 1). Since 
both p(r) and f(r) are real functions, the Fourier 
transform of ip(r)f(r) is purely imaginary. We 
have therefore demonstrated that in this special 
configuration, the aberration-free CLOWFS com- 
plex amplitude distribution in the CLOWFS de- 
tector array and the change introduced on this 
complex amplitude by a focus aberration are per- 
fectly orthogonal. This problem can be solved by 
introducing a defocus in the CLOWFS re-imaging 
optics, as proposed in this paper. This defocus 
term is equivalent to convolving the complex am- 
plitude map in the annulus by a kernel which 
breaks the orthogonality. 

Numerical simulations show that the CLOWFS 's 
response is linear to the modes considered in this 
paper. This linearity is only valid for small ampli- 
tudes, as shown in Figure [5] for tilt. 

4.2. Wavefront control algorithm and im- 
plementation 

As shown in Figure [HI the main steps of the 
CLOWFS closed loop are : 

• Acquire a single frame with the CLOWFS 
detector array 

• Compute the difference between this frame 
and the " reference" image obtained (or com- 
puted) with no aberrations. 
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Tip aberration (lambda/D) 

Fig. 5. — Linearity of the CLOWFS response to 
a tilt aberration. With a linear wavefront re- 
construction algorithm, the measured Tip signal 
(solid line) differs slightly from the true Tip aber- 
ration (dashed line). This non-linearity reaches 
about 20% at 0.2 A/D. 

• Decompose this difference as a linear sum 
of modal responses. These modal responses 
are either pre-computed or measured prior 
to starting the CLOWFS loop. 

• The coefficients of the linear decomposition 
described above are used to drive actuators 
to remove the aberrations measured by the 
CLOWFS 

The estimation of low-order aberrations there- 
fore requires knowledge of the modal responses 
(how aberrations will linearly modify the image 
acquired by the CLOWFS). As shown in Figure [6] 
(left), this information can consist of CLOWFS re- 
sponses to a series of wavefront aberration modes, 
which are here Zernike aberrations introduced ei- 
ther before (Ml to M5) or after (m6 to m8) the 
PIAA optics. They can either be pre-computed by 
simulation or measured as a response to an aber- 
ration introduced in the optical train (this latter 
method is more robust as it will accurately ac- 
count for unknown fabrication errors in CLOWFS 
optical components). 

The wavefront control algorithm must avoid 
confusion which could arise from the fact that 
the same low order aberration will produce nearly 
the same CLOWFS signal (although with a dif- 
ferent amplitude) wether it is introduced before 
or after the PIAA optics: in the configuration 
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Fig. 6. — Left: CLOWFS calibration prior to sensing. The CLOWFS response to 8 wavefront aberration 
modes (5 aberrations before the PIAA optics + 3 aberration after the PIAA optics) is measured or computed. 
The strong similarity between pre-PIAA and post-PIAA modes requires computing " differential aberration" 
responses M6, M7 and M8. The resulting basis of responses (Ml to M8) is orthogonal, and can then be 
used to decompose the CLOWFS image in a linear sum of aberrations (top right). Right: An example of 
a control scheme is shown for tip/tilt modes only (center/bottom right). In this Figure, both inside and 
outside focus CLOWFS images are shown. The response to Zernike aberrations for inside and outside focus 
are either equal (for antisymetric Zerike modes) or opposite (for symetric Zernike modes). The CLOWFS 
sensitivity is therefore identical whether the image acquired is inside or outside focus, or if both images are 
acquired. 



adopted in this paper, responses to pre-PIAA 
and post-PIAA tip/tilt are 99.6% similar. Al- 
though post-PIAA aberrations are expected to be 
smaller/slower, they still need to be properly cor- 
rected. One must avoid compensating for a tele- 
scope pointing error (pre-PIAA tip-tilt) by intro- 
ducing a post-PIAA tip-tilt, or compensating for a 
post-PIAA tip-tilt aberration by a pre-PIAA tip- 
tilt correction: both scenarios would create strong 
diffracted light at and beyond the IWA in the sci- 
ence focal plane even though the "overall" tip-tilt 
signal seen by the CLOWFS would be zero. The 
reconstruction algorithm shown in Figure [6] ad- 
dresses this issue by creating "differential" tip-tilt 
and focus modes (M6, M7, M8) obtained by sub- 
traction of pre-PIAA responses M1,M2, M3 from 
post-PIAA responses m6, m7, m8. These differ- 



ential modes track the difference between pre and 
post-PIAA aberrations, and can only be measured 
slowly since the corresponding CLOWFS response 
is weak (due to the high degree of similarity be- 
tween modes Ml, M2, M3 and m6, m7, m8). Since 
most tip-tilt /focus aberrations originate before the 
PIAA optics, the proposed algorithm offloads all of 
the tip/tilt signal to telescope pointing corrections 
(pre-PIAA). Differential aberrations are measured 
more slowly and off-loaded as post-PIAA correc- 
tions. In the example given here, thanks the large 
gap between pre-PIAA and post-PIAA pointing 
control bandwidths (w kHz vs. » Hz) the differ- 
ential pointing signal can be sent to a post-PIAA 
tip/tilt corrector rather than a combination of pre- 
PIAA and post-PIAA tip-tilt. 
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Fig. 7. — Tip (top left), focus (top right), astigmatism (bottom left) and differential tip (bottom right) 
sensitivity of the CLOWFS as a function of the relative size of the opaque disk in the focal plane mask 
(fiA"2) and the CLOWFS dcfocus distance. The sensitivity map is shown as a grey scale 2D map and the 
corresponding projection on the ri/r 2 and de focus axes are shown as plots above and to the left of each 2D 
map. Sensitivities are measured as the dispersion on a sample of 10 uncorrelated measurements with 10 6 
photons at the telescope entrance each, and are shown here scaled to one photon (equal to the dispersion 
multiplied by the square root of the number of photon). 
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4.3. Closed loop accuracy and stability 

The accuracy of the CLOWFS loop is ulti- 
mately limited by the measurement accuracy (as- 
suming ideal actuators). For small wavefront ex- 
cursions, the linear model is a very good approx- 
imation and the measurement accuracy is driven 
by detector noise and photon noise. For optimal 
performance, the CLOWFS exposure time and 
loop gain need to be carefully chosen to match 
the dynamic properties of the wavefront aberra- 
tions. In principle, in a very stable environment, 
the CLOWFS loop could correct very small er- 
rors by averaging a large number of measurements 
(small loop gain). Quantitative estimates will be 
given in Sj5] assuming a photon- noise limited detec- 
tor. 

The closed loop controller proposed in this pa- 
per is based on the linearity of the CLOWFS re- 
sponse. For large wavefront excursions, the actual 
CLOWFS response can differ significantly from 
this linear model, and the closed loop may "un- 
lock". While we have not performed a detailed 
loop stability analysis, we expect the CLOWFS 
loop to stay locked as long as pointing excursions 
are within a few tenths of X/D: 

• Figure [5] shows that the CLOWFS is still 
very close to the linear model at 0.2 X/D, 
and, more importantly, the curve shown in 
this figure is smooth and monotonic. With 
a 0.2 X/D pointing excursion, the estimate 
obtained from our linear model will not be 
exact, but will be sufficiently good to bring 
back the pointing significantly closer to on- 
axis at the next iteration 

• As illustrated in Figure HI the CLOWFS re- 
covers pointing errors by essentially looking 
for a brightness enhancement on one side 
of the re-imaged reflective ring and a cor- 
responding decrease in surface brightness on 
the opposite side of the ring. As long as 
pointing errors are small enough for the peak 
of the stellar image to stay within the in- 
ner edge of the reflective ring, this scheme 
will correctly estimate the direction of the 
pointing error, and the amplitude of this er- 
ror will be estimated with sufficient precision 
to reduce the pointing error after correction. 
Since the inner edge of the reflective mask 



is at least 0.5 X/D, the pointing loop is ex- 
pected to stay stable over this range. 

Laboratory operation of the CLOWFS confirm 
this behavior, and the closed loop locks even if 
the initial pointing error is 0.5 X/D. 

If the initial pointing error is too large for the 
linear CLOWFS loop to lock, a non-linear model 
of the CLOWFS may be used. Even with a non- 
linear model, the CLOWFS will fail to measure 
pointing errors if the stellar PSF misses the fo- 
cal plane mask (pointing error is as large or larger 
than the mask size). A separate coarse pointing 
sensor is then necessary to measure the pointing 
offset and then bring the stellar PSF within the 
CLOWFS range. The coronagraph image may be 
used for this purpose, provided that detector sat- 
uration does not prevent pointing error measure- 
ment. 

4.4. Optimal CLOWFS design for PIAA 
coronagraph 

In this section, we explore how the CLOWFS 
performance is affected by the choice of r\/r<i 
(which sets the amount of light sent to the 
CLOWFS) and the amount of de focus introduced 
in the image of the focal plane mask. For each 
choice of the parameters r\/ri and de focus, the 
linear algorithm described in S 14.2I is used here 
to measure the CLOWFS sensitivity to point- 
ing (modes Ml, M2), focus (mode M3), astigma- 
tism (modes M4, M5), differential tip-tilt (modes 
M6,M7) and differential focus (mode M8). The 
linear algorithm allows rapid evaluation of a large 
number of designs: 2-D sensitivity maps on a 
r\/ri - de focus plane can therefore be created, as 
shown in Figure [Tj 

Under the ideal conditions used in this simu- 
lation (perfect detector) , Figure [7] shows that all 
aberrations are best measured if r\/r<z is close to 
zero: ideally, the focal plane mask should be fully 
reflective. Figure[7]however shows that there is rel- 
atively little loss in sensitivity if r% / r-i is increased 
to up to w 0.4. With only 10% of the starlight 
reflected to the CLOWFS (n/r 2 = 0.294) there 
is almost no loss in CLOWFS sensitivity to point- 
ing errors. For the pointing error to double, T\/ri 
must be increased to 0.44, leaving only 0.87% of 
the starlight for the CLOWFS. If the same 0.87% 
were applied uniformly over the PSF (instead of 
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masking the center of the PSF), pointing errors 
would increase more than tenfold. A similar be- 
havior is observed for all other modes shown in 
Figure [7) for example, the focus measurement er- 
ror is doubled for ri/r 2 = 0.56, leaving only 0.06% 
of the starlight for the CLOWFS. These results 
confirm the "relative signal amplification" effect 
previously claimed: the central dark part of the 
mask removes most of the CLOWFS flux but only 
mildly reduces the signal. The CLOWFS design 
adopted in this paper (r 1 /r 2 — 0.4, de focus = 20 
rad PV) uses only 1.8% of starlight, at the expense 
of almost doubling tip-tilt and focus measurement 
errors. 

The amount of defocus has little effect on the 
CLOWFS sensitivity to pointing, unless it is suffi- 
ciently large (> 40 rad) to "blend" together light 
from opposite sides of the reflective focal plane 
ring. Some defocus is however essential to al- 
low good focus estimation with the CLOWFS, and 
helps with astigmatism as well. For focus, there 
is however no gain beyond defocus ~ 20 rad, at 
which point the loss in pointing sensitivity is still 
very small. 

For all pre-PIAA low order aberrations, Fig- 
ure [Jj shows that the residual sensing errors with 
a CLOWFS are within a factor «2 of the the- 
oretically optimal sensitivity. For example, the 
CLOWFS pointing error is « 1 rad per mode per 
photon, or w 1.4 rad for tip and tilt combined. 
This corresponds to just under 1 A/D pointing er- 
ror for a single photon, even if the central opaque 
zone of the focal plane mask removes most of the 
starlight. Similarly, Focus and astigmatisms are 
recovered with w 2 rad RMS error each for a single 
photon. The CLOWFS therefore makes a very ef- 
ficient use of a limited number of photons. Due to 
the strong similarity between pre-PIAA and post- 
PIAA modes, the CLOWFS sensitivity to differen- 
tial modes such as M6, M7 or M8, is much weaker. 

Beyond the global trends outlined above, Fig- 
ure [Jj also shows more subtle effects: diffraction 
effects produce oscillations of the sensitivity to Fo- 
cus and Astigmatism in the ri/ra - defocus plane. 
Computing 2D sensitivity maps such as the ones 
shown in FigureUJis therefore essential for fine tun- 
ing the CLOWFS performance to low order wave- 
front control requirements. 



5. Aberration sensitivity of a PIAA coro- 
nagraph equipped with a CLOWFS 

5.1. Comparison between CLOWFS and 
focal plane based wavefront control 

We consider here a 1.4-m diameter telescope 
observing a my = 6 star in a 0.2 /mi wide band 
centered at A — 0.55/im. We assume a 50% system 
throughput, which corresponds to 6 x 10 9 ph.s~ l 
at the telescope entrance. As described in Sj2] 
pre-PIAA low order aberrations can affect coro- 
nagraphic contrast at small angles (close to the 
IWA of the PIAA coronagraph) much more eas- 
ily than post-PIAA low order aberrations. We 
therefore only consider in this section pre-PIAA 
pointing errors, which could be generated by a 
telescope/spacecraft pointing error. We adopt the 
requirement defined in pointing error must be 
measured to a 0.0016 X/D 1-cr accuracy. 

Without CLOWFS, a 0.0016 A/D pointing er- 
ror would have to be measured from the corre- 
sponding 6.9 x 10~ n total coronagraphic leak in 
the science focal plane, equal to approximately 
0.4 ph.s -1 . A 1-cr measurement of this leak can 
be achieved with « 1 photon, provided that this 
leak is interfcromctrically combined with a much 
brighter and well known "reference(s)" (in 2-D, 
this technique is referred to a Focal plane wave- 
front sensing, where a Deformable Mirror is used 
to mix coherent starlight with the scattered light 
halo). This "reference(s)" is necessary to measure 
the sign of the aberration and also to bring the 
signal above the detector readout noise and/or in- 
coherent background in the image. In this scheme, 
detection of a faint companion cannot be done at 
the same time as pointing measurement, and time 
must be shared between the two tasks (we assume 
here that 50% of the time is spent for each task). 
A measurement of pointing error with a 0.0016 
X/D 1-cr error would therefore require 5 seconds. 

As shown in Figure [Jj the tip sensitivity for the 
CLOWFS is 1.2 rad RMS for one photon at the 
telescope entrance. The 0.0016 X/D tip measure- 
ment accuracy (equal to 0.0025 rad RMS) there- 
fore requires 2.3 x 10 5 photons, which can be gath- 
ered in 38 /is. The CLOWFS is therefore capable 
of measuring pointing errors about 130,000 times 
faster than the science focal plane. 

The pointing stability requirement is derived in 
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Table 1 

Pointing stability requirements for a PIAA coronagraph with and without CLOWFS 



Without CLOWFS With CLOWFS 



Required pointing calibration accuracy (10 11 contrast) 0.0016 A/D (0.13 mas) 

Maximum RMS pointing excursion (10 -10 contrast) 0.005 A/D (0.4 mas) 

Required sampling time b 5 s c 38 fj,s 

Maximum allowed uncalibrated pointing drift rate 0.026 mas.s' 1 3.4 arcsec/s 



a For detection of a 10 10 contrast source around a my = 6 star observed with a 1.4-m telescope in a 0.2/xm wide 
band centered at 0.55 /im with a 50% system throughput. 

b Sampling time required to measure the pointing error with a 1-a error equal to the "Required pointing calibration 
accuracy" . 

c Assumes that 50% of the observing time is dedicated to measurement of low order aberrations. Also assumes that 
the signal is well above readout noise and zodi/exozodi background. 



table [T] from both the required pointing calibra- 
tion accuracy (equal to 0.0016 X/D in this exam- 
ple) and the sampling time necessary to measure 
pointing errors to this level of accuracy. The ratio 
between these two quantities defines a maximum 
allowable pointing drift rate beyond which the sta- 
bility/calibration requirement previously defined 
cannot be met. As shown in Table [TJ measuring 
pointing errors with the CLOWFS is several orders 
of magnitude quicker than if only the light in the 
science focal plane were used. This measurement 
sensitivity gain yields a much more relaxed point- 
ing drift rate requirement: SAarcsec.s -1 with a 
CLOWFS as opposed to 0.026 mas. s^ 1 without. 

5.2. Detector requirements 

The CLOWFS camera requires a modest num- 
ber of pixels: the signals shown in the left of Figure 
[6] contain little high spatial frequencies and can be 
accurately measured with approximately 10 pix- 
els across one defocused spot images. A 20 by 20 
pixel window (400 pixels) is sufficient, and can be 
read rapidly (>10 kHz) with current technology. 

The CLOWFS measures pointing and other low 
order aberration by detecting changes in the defo- 
cused image it acquires. Temporal changes in the 
detector response must therefore be small com- 
pared to the expected signal. The CLOWFS is not 
sensitive to static spatial variations in the detector 
response (flat field) , as the signal is extracted from 
a difference between two images acquired at differ- 



ent times. The effect of spatially uniform flat field 
variations can be removed by scaling of the images 
prior to this subtraction. A 0.0016 A/D pointing 
offset corresponds to a 2% change in the surface 
brightness on the bright ring images by each frame 
of the CLOWFS: one side of the ring is 2% brighter 
while the opposite side is 2% fainter. A compara- 
ble variation in detector sensitivity between the 
two sides of the ring images is very unlikely in 
modern detectors (visible CCDs), even over the 
course of several hours. 

Thanks to the large number of photons col- 
lected by the CLOWFS, a moderate amount of 
readout noise (few photo-electron) will still allow 
operation of the CLOWFS at high sampling rate. 
For example, at 10 kHz sampling rate, the defo- 
cused image contains 10800 photons on a my = 6 
target. Assuming that the "ring" occupies a 50 
pixel surface area on the detector, photon noise 
(15 photo-electrons per pixel) is likely to be larger 
than readout noise on modern visible detectors. 

6. Laboratory Demonstration 

The CLOWFS has been implemented in the 
PIAA coronagraph testbed at the Subaru Tele- 
scope. As shown in Fig. the CLOWFS is driv- 
ing 5 piezo actuators to move the light source in 
x,y,z and move a post-PIAA mirror in tip and tilt. 

The CLOWFS focal plane mask used in this 
experiment is shown in Fig. [9l and was manu- 
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Fig. 8. — Schematic representation of the 
CLOWFS implementation for the PIAA coro- 
nagraph testbed at Subaru Telescope. The 
CLOWFS is extracting light reflected by the fo- 
cal plane mask and drives 5 actuators: pre-PIAA 
tip/tilt /focus (the light source is mounted on a 3 
axis piezo stage) and post-PIAA tip/tilt. 

factured by lithography techniques. The mask is 
not "ideal": the central portion is not perfectly 
opaque, and reflects a few percent of the light. 
Since most of the starlight falls on this central 
part of the mask, the CLOWFS frame contains 
a bright peak at the center of the defocused mask 
image. Thanks to the defocus introduced in the 
CLOWFS, the fainter light reflected by the reflec- 
tive annulus interferes with this central peak, and 
also forms fainter outer rings visible as shown in 
Fig. [TO] With a finite reflectivity of the central 
part of the focal plane mask, a more optimal de- 
sign would be to reduce the size of the central 
"dark" area to increase the amount of light at the 
transition between the central dark spot and the 
reflective annulus. 

The linear scheme implemented for calibration 
and processing of the CLOWFS data is very in- 
sensitive to the details of the CLOWFS design, 
and the CLOWFS is able to measure simultane- 
ously both pre and post-PIAA tip/tilt with little 
cross-talk (see Fig. [TUl upper right). As shown in 
Fig. EH we have achieved ~ 10 _3 A/D closed loop 
pointing stability for both the light source position 
and post-PIAA tip/tilt. 

In our laboratory testbed, pre-PIAA pointing 
errors and post-PIAA tip/tilt are similar in am- 
plitude, while in an actual system, we would ex- 
pect pre-PIAA pointing errors (telescope pointing, 
primary and secondary mirrors tilts) to be much 
larger and faster than post-PIAA tip/tilt (most 
likely due to slow thermal drifts on small optics). 



Fig. 9. — CLOWFS focal plane mask used in the 
PIAA coronagraph laboratory testbed at Subaru 
Telescope (fabricated by HTA photomask). The 
100 micron radius mask center is opaque (low re- 
flectivity), and is surrounded by a 100 micron wide 
highly reflective annulus. The science field, trans- 
miting light to the science camera, extends from 
200 micron to 550 micron radius. 

On our testbed, we have therefore operated the 
control loop with similar temporal bandwidth for 
both pre and post-PIAA aberrations, unlike the 
control scheme proposed in Fig. [6l 

The performance we have achieved in the lab- 
oratory is limited by both system stability (our 
testbed is in air and includes 75 mm diame- 
ter optics separated by more than a meter) and 
CLOWFS loop speed. In our experiment, the 
CLOWFS sampling interval was limited to 25 s 
due to hardware limitations (readout speed of 
the camera image and time necessary to trans- 
fer the image from the computer which controls 
the camera to a separate computer performing the 
CLOWFS image analysis). For a space corona- 
graphic mission, a better controlled environment 
and a faster readout camera (10 kHz is reasonable 
for the small number of pixels needed) would allow 
higher performance. Despite these limitations, our 
laboratory demonstration of the CLOWFS con- 
cept has exceeded both the RMS pointing excur- 
sion and the pointing measurement accuracy re- 
quired for achieving lelO contrast at 2 X/D with 
a visible PIAA coronagraph mission. 
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Fig. 10. — Laboratory performance for the CLOWFS. Upper left: Measured CLOWFS reference frame 
and influence functions for the 5 axis controlled in the experiment. Pre-PIAA and post-PIAA modes look 
extremely similar, as expected. Top right: Open loop simultaneous measurement of pre and post-PIAA 
modes. The measured amplitudes match very well the sine-wave signals sent to the actuators, and the 
CLOWFS is able to accurately measure all 4 modes shown here with little cross-talk. Since this measurement 
was performed in open loop, the measurement also include unknown drifts due to the limited stability of 
the testbed. Bottom left: Closed loop measurement of the residual error for the 5 modes controlled. The 
achieved pointing stability is about 10~ 3 X/D for both the pre-PIAA and post-PIAA tip/tilt. Bottom right: 
Position of the actuators during the same closed loop test. 



7. Conclusion 

The CLOWFS design presented in this paper 
can efficiently measure low order aberrations "for 
free" , as it uses light that would otherwise be dis- 
carded by the coronagraph. Both the hardware 
configuration and software algorithms presented 
are easy to implement and their performance is ro- 
bust against calibration errors, chromaticity, non- 
common path errors and small errors/aberrations 
in the optical components. The CLOWFS point- 
ing measurement can also lead to improved as- 
trometric accuracy for the position of faint com- 
panions, as the star position on the image is usu- 



ally difficult to me asure in coronagraphic images 
(|Digbv et alJl2006h - 



In this paper, we have studied a CLOWFS de- 
sign on a low-IWA PIAA coronagraph. Although 
coronagraphs with larger IWA can tolerate larger 
amount of low order aberrations (see for example 
Kuchner & Traub 2002; Kuchner et al. 2005), they 
require these aberration to be measured ahead of 
the coronagraphed beam because they are " blind" 
to aberrations until they are large enough to pro- 
duce significant coronagraph leaks. The CLOWFS 
would therefore be very useful to any high contrast 
coronagraph, and the optical design presented in 
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this paper can readily be used on any corona- 
graph where a focal plane mask physically blocks 
starlight. 

The CLOWFS can also be applied to phase 
mask coronagraphs. For such coronagraphs (see 
for example Roddier & Roddier 1997; Rouan et al. 
2000; Palacios 2005), where starlight is diffracted 
outside the pupil by the focal plane mask, a modi- 
fied Lyot stop is placed on the pupil plane to reflect 
starlight to the CLOWFS. In addition, using a 
pattern matching algorithm, the CLOWFS can es- 
timate low-order wavefront aberrations accurately 
and quickly even in non-linear region. These new 
CLOWFS designs and their performances will be 
presented in an upcoming paper. 
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